INVARIANTS ASSOCIATED WITH LINEAR ORDINARY 
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Abstract. We apply a novel method for the equivalence group and its infin- 
itesimal generators to the investigation of invariants of linear ordinary differ- 
ential equations. First, a comparative study of this method is illustrated by 
an example. Next, the method is used to obtain the invariants of low order 
linear ordinary differential equations, and the structure invariance group for 
an arbitrary order of these equations. Other properties of these equations are 
also discussed, including the exact number of their invariants. 



1. Introduction 

Due to promising results obtained on invariants of algebraic functions, and to 
the similarities in properties between differential equations and algebraic equations, 
the study of invariants of differential equations began in the middle of the nineteen 
century. One of the most important studies of these functions was carried out by 
Forsyth in his very valuable memoir [5], in which he considers the linear ordinary 
differential equation of general order n in various canonical forms, the first of which 
is given by 

+ a„_iy(""i) + a„-22/^"-'' + ■ ■ ■ + aoy = 0, (1.1) 

where the coefficients a„_i, a„_2, . . . , ao ai'e arbitrary functions of the independent 
variable x. Earlier writers on the subject, cited here in an almost chronological or- 
der, include Laplace [10 , Laguerre "9], Brioschi 1], and more importantly Halphen, 
who made ground-breaking contributions in his celebrated memoir 3 . 

Methods used up to the middle of the nineteen century for studying invariant 
functions were very intuitive, and based on a direct analysis, in which most results 
were obtained by a comparison of coefficients of the equation before and after it was 
subjected to allowed transformations. However, the application of these techniques 
has been restricted to linear equations. 

Based on ideas outlined by Lie [llj . the development of infinitesimal meth- 
ods for the investigation of invariants of differential equations started probably 
in [18] ■ and a formal method has been suggested 4, which is being commonly 
used [m m [6l m [7] . However, the latter method requires the knowledge of the 
equivalence transformations of the equation, and thus a new method that provides 
these transformations and at the same time the infinitesimal generators for the 
invariant functions has recently been suggested [14]. 
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In this paper, we use the method of [1^ to derive exphcit expressions for the 
invariants of various canonical forms of the general linear ordinary differential equa- 
tions of order up to 5. We also use the same method to derive the structure invari- 
ance group for a certain canonical form of the equation, and for a general order. 
Relationships between the invariants found as well as some of their other proper- 
ties, including their exact number, are also investigated. We start our discussion 
in Section [2] by an illustrative example comparing the former infinitesimal method 
of nig with that of [H]. 

2. Two METHODS OF DETERMINATION 

We begin this section with some generalities about equivalence transformations. 
Suppose that T represents a family of differential equations of the form 

A(x,2/(„);C) = 0, (2.1) 

where x = {x^, . . . , x") is the set of independent variables, y = {yi, . . . , yq) is the 
set of dependent variables and j/(„) denotes the set of all derivatives of y up to the 
order n, and where C denotes collectively the set of all parameters specifying the 
family element in T. These parameters might be either arbitrary functions oi x, y 
and the derivatives of y up to the order n, or arbitrary constants. Denote by G a 
connected Lie group of point transformations of the form 

X ^ (j){x,y](j), y ^ip{x,y;(T), (2.2) 

where a denotes collectively the parameters specifying the group element in G. 
We shall say that G is the equivalence group of (|2.1|) if it is the largest group of 
transformations that maps every element of J- into In this case (|2.2|) is called the 
structure invariance group of (|2.ip and the transformed equation takes the same 
form 

A(x,y(„);C) =0, (2.3a) 

where 

C, =Q(C,C(,);a), (2.3b) 

and where C(s) represents the set of all derivatives of C up to a certain order s. In 
fact, the latter equality (j2.3b[) defines another group of transformations Gc on the 
set of all parameters C of the differential equation [12], and we shall be interested 
in this paper in the invariants of Gc- These are functions of the coefficients of the 
original equation which have exactly the same expression when they are also formed 
for the transformed equation. We also note that G and Gc represent the same set, 
and we shall use the notation Gc only when there is need to specify the group action 
on coefficients of the equation. 

The terminology used for invariants of differential equations and their variants 
in the current literature [13 [6l [H [7l [13] is slightly different from that of Forsyth [2] 
and earlier writers on the subject. What is now commonly called semi-invariants 
are functions of the form $(C, C(r)) that satisfies a relation of the form $(C, C(r)) — 
w(cr) • $((7, C'(^)). When the weight w(cr) is equal to one, <i>(C, C(j.)), is called an 
invariant, or an absolute invariant. Semi-invariants usually correspond to partial 
structure invariance groups obtained by letting either the depend variable or the 
independent variables unchanged. 

Because the infinitesimal method proposed in [14] is still entirely new and has 
been, in particular, applied only to a couple of examples, we wish to illustrate 
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a comparison of this method with the most well-known method of For this 
purpose we consider Eq. with n = 3, which is the lowest order for which a 

linear ODE may have nontrivial invariants. The structure invariance group of p.ip 
is given by the change of variables x = f{x), y = T{x)y, where / and T are 
arbitrary functions, and for n ~ 3, this equation takes the form 

y^^^ + a^v" + aiy' + a^y - 0, (2.4) 

and it is easy to see that the corresponding infinitesimal transformations associated 
with the structure invariance group may be written in the form 

X K. X -^r ei{x), y K y + e{i]{x)y), (2.5) 

where ^ and rj are some arbitrary functions. If we now let 

K^'^ ^^d.,+yVdy+Cldy, +C2 dy„ + Cs Oy.. 

denote the third prolongation oi Va — dx + yrj dy, then we may write 



y' « y' + eCi{x,y,y'), y" ^ y" + e C2{x,y,y' ,y"), y « y'^' + £ (3(2;, y(3))- 

/2.6) 

According to the former method [Ij, in order to find the corresponding infin- 
itesimal transformations of the coefficients a2 , ai , and oq of (|2.4p , and hence to 
obtain the infinitesimal generator for the associated induced group Gc, the original 
dependent variable y and its derivatives are to be expressed in terms of y and its 
derivatives according to approximations of the form 

y' Ky' -eCi{x,y,y'), y" ^ y" - e C2{x,y,y' ,y"), y'^) k y'^^) - e C-,{,x,y,y(^s)). 

(2.7) 

These expressions are then used to substitute y and its derivatives for y and its 
derivatives in the original differential equation. In certain cases such as that of 
simple infinitesimal transformations of the form (j2.5p . the required approxima- 
tions (|2.7p can be obtained from (|2.6[) . More precisely, an explicit calculation of 
Va^^ shows that 

Ci = yg' + (.9 - />' (2.8a) 
C2 - y'(2.g' - /") + yg" + [g - 2/')y" (2.8b) 

C3 = 3(5' - nv" + y'iig" - n + yg'" + ig- 'W (2.8c) 

The first explicit approximation is readily obtained from the second equation 
of (j2.5|) which shows that 

y = y-^^y{l-eg), (2.9) 
1-1- £5 

by neglecting terms of order two or higher in e. In a similar way, using the equa- 
tions (|2.7p and (|2.8p we obtain after simplification the following approximations 



V' = {y' ~eyg'){l~e{g- f')) (2.10a) 
y" ^[y"~e {y'{2g' - /") + yg") (1 - s{g - 2/'))] (2.10b) 

y'" = [y - e (3(g' - /")y" + y\ig" - /"') + yg^^^)] (l - e{g - 3/')) (2.10c) 
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Substituting equations (12. 9|) and ()2.10p into ()2.4p and rearranging shows that 
the corresponding infinitesimal transformations for the coefficients of the equation 
are given by 

do = ao + £[— 3ao/' - aig' - a2g" - g'^'] (2.11a) 

di = fli +e[-2aif' - 2a2.g' + fla/" - 3.g" + J^^)] (2.11b) 

a2 =a2 + e[-a2/'-35' + 3/"]. (2.11c) 

In other words, the infinitesimal generator of the group Gc corresponding to (|2.4p 
is given by 

+ (-02/' - 3.9' + 3/") da, 

(2 12) 

+ (-2ai/' - 2a25' + fl2/" - 3.g" + /(3)t) 

+ (-3ao/' - ai.9' - a2.g" - .9^'') 9aa 

The other method that has just been proposed in |14j for finding the infinitesi- 
mal generators of the group Gc for a given differential equation simultaneously 
determines the equivalence group of the equation, in infinitesimal form. The first 
step in the determination of X'^ is, according to that method, to look for the infini- 
tesimal generator X of ()2.4|) , in which the arbitrary coefficients are also considered 
as dependent variables. For a general equation of the form (|2.ip . X takes the form 



X = ^idxi ^ ^^pdxp +'nidyi^ \-Vq dy^ + (t>i del + 02 <9c2 + . . . , 

where (C^, C^, . . . ) = C denotes the set of parameters specifying the family element 
in J-'. This expression for X may be written in the shorter form X = , 77, 0} , 
where ^,ry and 4> represents collectively the functions and (/>j, respectively. 

The next step according to the method is to look at minimum conditions that 
reduce V — {C,??} to a generator V'^ — {C°7^°} of the group G of equivalence 
transformations. These conditions are also imposed on </> to obtain the resulting 
function (fP, and hence the generator = {^°, (/i^} of Gc. 

In the actual case of equation (|2.4p we readily find that the generator 

X = ^dj; + ?]dy + (j)2 da, + 01 dai + 00 dao (2.13) 

is given by 

i = f 

v = 9y + h 

02 -02/' - 35' + 3/" 

01 = -2ai/' - 2a25' + ^2/" - 3^" + 

00 = aoh + iaoyf + aiyg' + aih'x + a2yg" + a2h" + yg'^'^^ +h^^^ , 

y I ■ J 

where /, g and h are arbitrary functions of x. It is clear that, due to the homogeneity 
of (|2.4p . we must have h = 0. This simple condition immediately reduces V = 
{/, yg + h} to the well-known generator V'^ = {/, yg} of G. If we denote by 0° = 
{02,0]',0o} the resulting vector when the same condition /i = is applied to 
= {02, 01, 00 } , then = {/, 0°} , which may be represented as 

X" = /a, + 0°9a.+0;aa, +0O9,„, (2.14) 
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is the same as the operator X'^ of (I2.12|) . 

It should however be noted that the former method involves a great deal of 
algebraic manipulations, and in the actual case of Eq. (|2.4p . the method works 
because of the simplicity of the infinitesimal generator Va = {/, yg} ■ For higher 
orders of the Eq. (jl.ip . difficulties with this method become quite serious. In fact, 
variants of this former method are being increasingly used, which also exploit the 
symmetry properties of the equation but even these variants still require the 
knowledge of the equivalence group. 



We shall use the indicated method of 14J in this section to derive explicit ex- 
pressions for invariants of linear ODEs of order not higher than the fifth. It should 
be noted that, not only the infinitesimal method we are using for finding these in- 
variants is new, but also almost all of the invariants found are appearing in explicit 
form for the first time. The focus of Forsyth and his contemporaries was not only 
on linear equations of the form (|2.ip . but also on functions $((7, C(r)) satisfying 
conditions of the form 



= or <D(C,CM) = (d5/da;r (3.1) 



in the case of linear ODEs, where h is an arbitrary function, cr is a scalar and where 
g — g{x) and x are defined by a change of variables of the form x = f{x), y = gy, 
where / is arbitrary. Such functions are clearly invariants of the linear equation only 
if the dependent variable alone, or the independent variable alone, is transformed, 
but not both as in our analysis. However, in his very talented analysis starting 
with the investigation of semi- invariants satisfying the second condition of (|3.ip . 
Forsyth obtained in his memoir [2] a very implicit expression for true invariants of 
linear ODEs of a general order, given as an indefinite sequence in terms of semi- 
invariants. As it is stated in [2], earlier works on the subject had given rise only 
to semi-invariants of order not higher than the fourth, with the exception of two 
special functions which satisfy the condition of invariance for equations of all orders. 

3.1. Equations in normal reduced form. Invariants of differential equations 
generally have a prominent size involving several terms and factors, and thus they 
are usually studied by first putting the equation in a form in which a number of 
coefficients vanish. By a change of variable of the form 



3. Invariants for equations of lowest orders 




(3.2) 



equation (jl.ip can be reduced to the simpler form 



y 



(3.3) 



after the renaming of variables, and where the Aj — Aj{x) are the new coefficients. 
Eq. (|3.3p . which is referred to as the normal form of (jl.ip . will be as in [2] our 
starting canonical form for the determination of invariants. 
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The generic infinitesimal generator X'^ of Gc suclr as tire one found in ()2.12|) 
linearly depends on free parameters Kj, and can be written in terms of these pa- 
rameters as a linear combination of the form, 



X"^J2^i^^' (3.4) 



where the Wj are much simplified vector fields free of arbitrary parameters and with 
the same number of components as X'^, and a function F — F{C) satisfies the condi- 
tion of invariance X^ • F = if and only if Wj ■ F — 0, ioi j ~ 1, . . . , i'. Consequently, 
the invariant functions are completely specified by the matrix M = {Wi, . . . , W^} 
whose jth row is represented by the components of Wj , plus the coordinate system 
in which the vector fields Wj are expressed. These considerations also apply to pro- 
longations of X° and corresponding differential invariants. For a given canonical 
form of the equation, we will use the symbols ^'J^''^ and X™ to represent explicit ex- 
pressions of invariants and corresponding infinitesimal generator respectively. 
In such a representation, the subscript a will denote the order of the equation, while 
the superscript b will represent the order of prolongation of the original generator 
X^, and c will represent the actual number of the invariant in some chosen order. 
The superscript m will represent the canonical form considered, and will be n in 
the case of the normal form (|3.3p . s for the standard form (jl.ip . and w for another 
canonical form to be introduced below. For consistency, coefficients in all canonical 
forms considered will be represented as in (|l.ip by the symbols Uj ~ o,jix)- 

For n = 3, the first nontrivial differential invariants of (|3.3p occurs only as from 
the third prolongation of This third prolongation of the generator X° has ex- 
actly one invariant V^s'^, and both X'^ = X^ and Vl/g'^ have already appeared in the 
recent literature [14]. However, we recall these results here in a slightly simplified 
form, together with those we have now obtained for the fourth prolongation of X^, 
as well as for the case n — A. 

For n = 3, we have 



X3" = fd^~ 2(ai/' + f^^^)da, + (-3ao/' - ai/" - f^^^) da, (3.5a) 

3,1 _ (9aiAi2 + 7M^-6/XM")' 
1728Ai8 



,T,3.1 \au,i^ -r I fJ. — u/^p ; /o ri \ 



= (3.5c) 

4.2 



*3 



18// 



(216a;5 - 324ag a[ + 18 ag(9af + 2ai^') + +9//^^^^ 



IS/x^ 



(/i'(28Ai' 2 + 9a; (ai a[ - 4/x")) - 9ao(3af -I- 4ai a>' ~ Sfi'fi"),) 



(3.5d) 



where / is an arbitrary function and where /i — —2ao + a^. For n = 4, the first 
nontrivial invariant occurs only as from the second prolongation of X'^ = X^, and 
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we have 

^4" = / a. + (-202/' - 5/(3)) a,, + (-3ai/' - 202/" - 5/(4)) 



^ -8ao/'-3(ai/" + a2/^'^) 



2 

= -(-lOOao + Qal + 50a'i - 20^') (3.6b) 

*4 ^ = Aii(-1053a^ - 4500a2(2a? - ai^ - 4^) + 180a^(130ao - 3(5a'i + 8 4')) 

- 100(1300a^ + 75fi{~10a'o + 3a") + 90a[a2 + 27a'p - 60ao(5a'i + 8 4')), 

(3.6c) 

where = ai - 4, and ^1 = 1/ [75000/^^^/^)] . 

3.2. Equations with vanishing coefficients a„_i and a„_2- A much simpler 
expression for the invariants is obtained if the two terms involving the coefficients 
a„_i and a,i_2 in Eq. (jl.ip are reduced to zero in the transformed equation. Such 
a change of variables can be accomplished by a transformation of the form 

12 f 

{x,x}^— — — — a„_2, y = exp(- / a„_idx)y, (3.7a) 

n[n-l){n+l) J 

where 

{x,x} = (x'x(3) - (3/2)x"2) x'-' (3.7b) 

is the Schwarzian derivative, and where x ' — dx/dx. Thus by an application of (|3.7p 
to (jl.ip . we obtain after the renaming of variables and coefficients an equation of 
the form 

2/^") + an-3J/'""'^ + «„-42/("-"^ + • • ■ + aoy = 0. (3.8) 

In fact, p.8p is the canonical form that Forsyth [2], Brioschi [T], and some of their 
contemporaries adopted for the investigation of invariant functions of linear ODEs. 
However, Forsyth who studied these equations for a general order did not derive 
any explicit expressions for the invariants, with the exception of a couple of semi- 
invariants. There is a number of important facts that occur in the determination 
of the invariants when the equation is put into the reduced form (13. 8p . but we 
postpone the discussion on the properties of invariants to the next Section, where 
a formal result regarding their exact number is also given. 

For n — 3, nontrivial invariants exist only as from the second prolongation of 
X° = X^, and we have computed them for the second and the third prolongation. 
For n = 4, there are no zeroth order invariants and we have computed these invari- 
ants for the first and the second prolongation of X™. For n = 5, the invariants are 
given for all orders of prolongation of X™, from to 2. Since invariants correspond- 
ing to a given generator are also invariants for any prolongation of this generator, 
its is only necessary to list the invariants for the highest order of prolongation of 
any given generator. In the canonical form (|3.8p . and regardless of the order of the 
equation, the generator of Gc will depend on three arbitrary constants that we 
shall denote by ki, k2 and ^3. 



9ao 



(3.6a) 
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For n = 3, we have 

3.2 _ i28a'i~36a„a'oa'^ + 9ala^^^) 

' - H 

For n ~ A, we have 

= [ki + x{k2 + k^x)] dx - 3ai(fc2 + 2kzx) da^ + [-3aifc3 - 4ao(fc2 + Sfcax)] d, 



ao 



IS _ _ (-2ao + a'l) 



"1 

/ / ao(ao+3a'i) n,3 



^2,2 ^ ^1,2 ^ v::!^ ^ 

2.3 _ (14a^ - Map + 3ai a'/)^ 
27af 

2 4 _ 16ag + 48a^ a!^ + %a\ a!^ - 9ao ai(6a[, + a![) 

FinaUy, for n — h we have 

X™ = [fci + .x(fc2 + k^x)\ dx - 3a2(fc2 + 2fc3a:) 3^, 

- 2 [-3a2fc3 + 2ai(fc2 + 2^30;)] da^ + [-4aiA;3 - 5ao(fc2 + 2k-ix)] dag 



(-ai + a^ 



M3 



2 3 3 (6a2a'j^ — a^ — 6ai aQ'^ 



*5' = *5 



216 a» 



^2,4 ^ ^1,4 ^ 5af + 9a^ a'p ~ 3ai a2(5ao + 2a'i) + 3af a^ 
^ ^ 9a^ ^ ■ 



and 

^2,5 ^ (7af - 14ai a'^ + 6a2 a^Q^ 

^ 216 a| 

^2 6 _ 4af + 24af 03 + 9a2a" - 9aia2(3ai + a'^ ) 

' ' 9^1 

2_7 [l8(-al - af + a^ a^') - 6aia^(lla;, + 2a'/) + af a2(55ao + 40ai + 6a^')] ^ 
^ ~ 5832 af 

3.3. Equations in the standard form (|l.ip . No attempt has ever been made 
to our knowledge, to obtain the invariant functions for equations in the canonical 
form (jl.ip . due simply to difficulties associated with such a determination and to 
the prominence in size that semi-invariants for this canonical form already display. 
Indeed, semi-invariants for the canonical form (jl.ip were calculated by Laguerre [3] 
for equations of the third order. Transformations of the form p.2p or p.7p are 
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usually applied to transform to an equation in which one or both of the coeffi- 
cients a„_i and a„_2 vanish. Even with the infinitesimal generator X° of ()2.12|) at 
our disposal, it is still difficult to find directly the corresponding invariants for the 
third order equation, which is the lowest for which the linear equation may have 
nontrivial invariants. However we show here by an example how these invariants 
can be found for the third order from those of equations in a much simpler canonical 
form. 

Under the change of variables p.2p . Eq. (12. 4p takes, after elimination of a con- 
stant factor, the form 

y(3)+Si^' + Soy = 0, (3.10) 

^^'^^^ Bo = (27ao-9aia2 + 2a^-94')/27 (3.11) 

Bi = (27ai - 9a^ - 274)/27. (3.12) 

The sole invariant ^'g'^ for equations of the form p.lOp and corresponding to the 
third prolongation of the associated generator X° is given by (|3.5p . If in that 
function we replace aq and oi by Bq and Bi respectively, and then express Bq and 
Bi in terms of the original coefficients ao,ai and 02, the resulting function is an 
invariant of (|2.4p . However, it does not correspond to the third, but to the fourth 
prolongation of the corresponding generator X° ~ that was obtained in (|2.12p . 
In other words, in the canonical form (|l.ip with n = 3, we have 

XI ^X", as given by ([TT2l) 

,T,4,i _ {9B^^^^ + 7^l'^-6^i^iy 

" 1728^1^ ' 

where we have fi = —2Bq + B'^. Using Theorem 2 of [14l, it can be shown that 
there are no invariants corresponding to prolongations of X| of order lower than 
the fourth, and that the fourth prolongation has precisely one invariant. 

4. A NOTE ABOUT THE INFINITESIMAL GENERATORS AND THE INVARIANTS 

The method of [TJ that we have used in the previous section provides the in- 
finitesimal generator of the equivalence group G for a given equation of a specific 
order. This infinitesimal generator must be integrated in order to obtain the corre- 
sponding structure invariance group for the given family of equations of a specific 
order. However, once the structure invariance group for such a given equation has 
been found, it is generally not hard to extend the result to a general order. We 
show how this can be done for equations of the form p.Sp . 

With the notations already introduced in Section [2l if we denote by X the full 
symmetry generator of (13. 8p for a specific order and then obtain the infinitesimal 
generator of the equivalence group G, we find that 

= Vo.s = [k2 + xih + k^x)] + y{ki + 2kix) dy, 

V° = Vb,4 = [k2 + x{k^ + k^x)] + |(2fci + 2k3 + 9kix) dy, 

V° = Vo,5 = [k2 + x{k3 + k4,x)] + y{ki + ks + Ak^x) dy, 

where the A:j,for j = 1, . . . ,4 are arbitrary constants. Upon integration of these 
three vector fields, we find that for n — 3,4,5, the structure invariance group can 



for n 


= 3, 


(4.1a) 


for n 


= 4, 


(4.1b) 


for n 


= 5, 


(4.1c) 
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be written in the general form 



'^'"""'^ y= .,. , (4.2) 



6(1 + ax) ' " + a.x)("-i) ^ 

where a, 6, c and d are arbitrary constants. To see why (|4.2p holds for any order n of 
Eq. (|3.8p . we first notice that these changes of variables are of the much condensed 
form 

x = g{x), y^T{x)y. 

It thus follows from the properties of derivations that when ()4.2|) is applied to (13. 8|) . 
every term of order m in y will involve a term of order at most m in y upon 
transformation. Since the original equation (|3.8p is deprived of terms of orders 
n—1 and n — 2, to ensure that this same property will also holds in the transformed 
equation, we only need to verify that the transformation of the term of highest order, 
viz. y*-"-*, does not involve terms in y("~i' or y("~^). However, it is easy to see that 
under (1121), y(") is transformed into 6"(1 + ax)"-~^^y^"-' /d. We have thus obtained 
the following result. 

Theorem 1. The structure invariance group of linear ODEs in the canonical 
form p.Sp . that is, in the form 

y(") + a„-3y'""') + a„-42;("-*^ + • • ■ + aoy = 0, 
is given for all values of n > 2 by (|4.2p . 

It is easy to prove a similar result using the same method for equations in the 
other canonical forms Ijl.ip and (j3.3p , for which the corresponding structure invari- 
ance groups are well-known. It is also clear that (|4.2p does not violate in partic- 
ular the structure invariance group of (|3.3p which is given [19] by a; = F(x) and 
y — aF'^'^~^'>^'^y, where F is an arbitrary function and a an arbitrary coefficient. 

All the invariants that we have found for linear ODEs are rational functions of 
the coefficients of the corresponding equation and their derivatives, and this agrees 
with earlier results j2] . Another fact about these invariants is that for equations in 
the form p.Sp . nontrivial invariants that depend exclusively on the coefficients of 
the equation and not on their derivatives, that is, invariants such as that in (|3.9ap 
which are zeroth order differential invariants of X'^ occurs as from the order 5 
onwards. Indeed, since X'^ depends only on three arbitrary constants in the case 
of the canonical form (|3.8p . there will be more invariants, as compared to the other 
canonical forms in which X° depends on arbitrary functions. More precisely, we 
have the following result 



Theorem 2. For equations in the canonical form (13. 8p . the exact number n of 
fundamental invariants for a given prolongation of order p > of the infinitesimal 
generator A"" is given by 

p(7i-2)-4, z/(n,p)^(3,0) ^^g^ 
" otherwise. 



Proof. We first note that the lowest order for equations of the form (|3.8p is 3, and 
we have seen that when n = 3,4, 5 the generator X'^ depends on three arbitrary 
constants alone, that we have denoted by ki, k2 and fca, and do not depend on any 
arbitrary functions. Consequently, this is also true for any prolongation of order 
p of X°, and it is also not hard to see that this property does not depend on the 



11 



order of the equation. It is clear that the corresponding matrix A4 whose jth row 
is represented by the components of the vector field Wj of ()3.4p will always have 
three rows and maximum rank r = 3 when {n,p) ^ (3, 0). On the other hand, X'"^ is 
expressed in a coordinate system of the form \x, a„_3, a„_4, . . . , ao} that has n — 1 
variables. Its pth prolongation is therefore expressed in terms of il/ — n— l+p(n— 2) 
variables. For ^ (3, 0), the exact number of fundamental invariants is M — r, 

which is n + "pin — 2) — 4, by Theorem 2 of |l4j. For = (3,0) a direct 

calculation shows that the number of invariants is zero. This completes the proof 
of the Theorem. □ 

Similar results can be obtained for equations in other canonical forms, but as it 
is customary to choose the simplest canonical form for the study of invariants, we 
shall not attempt to prove them here. 

It should be noted at this point that invariant functions have always been inti- 
mately associated with important properties of differential equations, and in partic- 
ular with their integrability. Laguerre [9] showed that for the third order equation 
(|2.4p . when the semi-invariant /i = —2Bo + B[ given in Eq. (|3.13p vanishes, there is 
a quadratic homogeneous relationship between any three integrals of this equation, 
and consequently its integrability is reduced to that of a second order equation. 
Subsequently, Brioschi [l] established a similar result for equations in the canonical 
form (|3.3|) . More precisely, he showed that for third order equations, when the 
semi-invariant /i = — 2ao + a'l of Eq. ()3.5|) vanishes, the equation can be reduced in 
order by one to the second order equation 

where y = and he also gave the counterpart of this result for fourth order 
equations when the corresponding semi-invariant fi — ai — a'2 oi Eq. (|3.6p vanishes. 
In fact all reductions of differential equations to integrable forms implemented in |3] 
are essentially based on invariants of differential equations. More recently, Schwarz 
[l9] attempted to obtain a classification of third order linear ODEs based solely on 
values of invariants of these equations. He also proposed [20| a solution algorithm for 
large families of Abel's equation, based on what is called in that paper a functional 
decomposition of the absolute invariant of the equation. 

5. Concluding Remarks 

We have clarified in this paper the algorithm of the novel method that has 
just been proposed in [14j for the equivalence group of a differential equation and 
its generators, and we have shown how its application to the determination of 
invariants differs from the former and well-known method of [4j [5] , by treating the 
case of the third order linear ODEs with the two methods. We have subsequently 
obtained some explicit expressions of the invariants for linear ODEs of orders three 
to five, and discussed various properties of these invariants and of the infinitesimal 
generator X". 

Another fact that has emerged concerning the infinitesimal generator of the 
group Gc for linear ODEs is that, whenever the equation is transformed into a 
form in which the number of coefficients in the equation is reduced by one, the 
number of arbitrary functions in the generator X*^ corresponding to the transformed 
equation is also reduced by one, regardless of the order of the equation. Thus while 
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depends on two arbitrary functions, X" depends only on one such function 
and X"^ depends on no such function, but on three arbitrary constants. This is 
clearly in agreement with the expectation that equations with less coefficients are 
easier to solve, since less arbitrary functions in X^ means more invariants and 
hence more possibility of reducing the equation to a simpler one. However, the full 
meaning of the degeneration of these functions with the reduction of the number 
of coefficients of the equation is still to be clarified. We also believe that on the 
base of recent progress on generating systems of invariants and invariant differential 
operators (see |16l[T7] and the references therein), it should be possible to treat the 
problem of determination on invariants of differential equations in a more unified 
way, regardless of the order of the equation, or the order of prolongation of the 
operator X'^ . 
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